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Abstract 



We use a generalised version of the individual-based Tangled Nature model of 
evolutionary ecology to study the relationship between ecosystem structure and 
evolutionary history. Our evolved model ecosystems typically exhibit interaction 
^ O , networks with exponential degree distributions and an inverse dependence between 

\ connectance and species richness. We use a simplified network evolution model to 

qh' demonstrate that the observed degree distributions can occur as a consequence of 

partial correlations in the inheritance process. Futher to this, in the limit of low 
connectance and maximal correlation, distributions of power law form, P{k)(xl/k, 
^ can be achieved. We also show that a hyperbolic relationship between connectance 

and species richness, C^l/D can arise as a consequence of probabilistic constraints 
, on the evolutionary search process. 

OO 
O 

\o 
o 

•2 : 1 Introduction 

The network topology of community structure is an important property of an ecosystem. 
. ^ ' plays an integral role in system dynamics, influencing future behaviour whilst emerging 

^ . as a consequence of dynamical history. This duality of cause and effect exemplifies the 
nonlinearities inherent to such systems and leads to differing theories about the origins 
of observed network properties. At least since May's work in the seventies (May, 1974) 
there have been attempts to explain the degree distributions and connectances of ecolog- 
ical networks by reference to stability issues (McCann, 2000) (Pimm, Lawton, & Cohen, 
1991) (Warren, 1994) (McCann, 2000). Here we focus on the the role played by correlated 
inheritance in the structuring of the interaction network and elucidate such effects by 
presenting two related models of evolutionary ecology. 

The first, the correlated Tangled Nature model, is an individual-based model of co- 
evolution in which the reproduction rate of the individual depends on interactions with 
co-existing types. Mutations move the population around in trait space, but an offspring 
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that differs by only a few mutations from the parent experiences a set of interactions very 
similar to that of the parent. The correlation in this model allows species to emerge as 
localised clusters of individuals in trait space. 

The second is a simplified version of the first and consists of a dynamical network with 
an invariant node number. At each time step a randomly selected node is removed from 
the network and replaced by a duplicate of another randomly selected node. During the 
duplication process the edges are inherited with a probability that ultimately determines 
the parent- offspring similarity. 

The correlated Tangled Nature model produces exponential degree distributions and 
we attribute this outcome to correlations in the inheritance process. We use the network 
model to demonstrate that correlated speciation dynamics can lead to such distributions 
and, additionally, to other forms observed in real ecological systems. Depending on the 
level of correlation, the network model distributions take binomial, exponential-like and 
power law forms with the latter tending to P{k) = 1/k in certain limits. Interestingly, 
the network model dynamics emulate ecological processes without any net growth in 
nodes nor explicit preferential attachment; aspects which are integral to many power- law 
network dynamics. 

The ensembled data from the Tangled Nature simulations also present another curi- 
ous network phenomenon: an inverse relationship between connectance and the number 
of species. This form of dependence has been observed in ecological field data; predomi- 
nantly food web structured but also mutualistic networks (Cohen & Briand, 1984)(01esen 
& Jordano, 2002)(Montoya & Sole, 2003)(Bascompte, Jordano, Melian, & Olcscn, 2003). 
In a generalised interaction network this relation is thought to represent a delineation of 
stable and unstable regions of network structure (May, 1974)(McKane, Alonso, & Sole, 
2000). Here, we formulate an alternative line of reasoning based on a balance of selection 
forces and the probabilistic limitations of phase space searching. A simple functional 
form is derived and then compared to data obtained from the simulations of the Tangled 
Nature model. 

In the following section we define the correlated Tangled Nature model and present 
the results regarding the degree distributions. To provide insight into these results and 
propose a causal mechanism for the distributional forms, we introduce the network evolu- 
tion model in section 4. Returning to the Tangled Nature model, in section 5 we present 
the second set of results concerning the inverse relationship between connectance and the 
number of species. A probabilistic model is then used to explain the phenomenon. 

2 Correlated Tangled Nature model 

This version of the Tangled Nature model succeeded its predecessor (Christensen, CoUo- 
biano. Hall, & Jensen, 2002) (Hall, Christensen, Collobiano, & Jensen, 2002)(Collobiano, 
Christensen, & Jensen, 2003) (Anderson & Jensen, 2005) by introducing correlations into 
the inheritance process. Here, we shall present the construct of the model whilst sum- 
marising the correlative effects. We refer the reader to the earlier paper for a more 
detailed explanation of how this is achieved (Laird & Jensen, 2005). 

Individuals, {a,P,...} are described by phenotype vectors of L = 16 traits, = 
(T", Tg-i.-.-i T^), with each trait taking a value from the periodically bounded range. 
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[0, 99999] . We then define a completely specified interaction matrix ( J-matrix) that 
accounts for all interactions between any possible phenotypes. Ultimately only a small 
number of distinct phenotypes will exist in our system and their interactions will be 
described by a small subset of the elements of this complete matrix. A proportion, 9, 
of the entries of the complete J-matrix, J(T",T'^), are assigned normally distributed 
values that have local correlations within the matrix structure. All other values of the 
remaining proportion of the matrix, 1 — 6, arc assigned zeroes which are treated as a lack 
of interaction between the two relevant phenotypes. The result is that given a mutation 
of one trait value we have an exponential decay in the correlation between parent and 
offspring interaction set values that is dependent upon the distance mutated in the trait 
value, 

C(J(T", T^), J(T/^, T^)) = ea;p[-A(T", T^)/^] e (0, 1], (1) 

Here, is the correlation length and C(J(T", T'^), J(T*^, T^)) is the correlation be- 
tween the interaction strengths of two phenotypes a and (3 when each are interacting with 
a third, 7. This measure shall be abbreviated to C(T", T'^) hereafter. If (3 is the mutated 
offspring of a, then it inherits the same phenotype vector but with a small random shift 
in one of the trait values. The function, A(T", T^) = | T^^^ — T^^J represents this single 
trait shift (when two similar phenotypes with multiple trait differences are considered, 
this function becomes more complicated; we refer the reader to the original paper for 
details (Laird & Jensen, 2005)). Importantly, interacting uncorrelated phenotypes take 
values J(T",T'') and J(T^, T") that are uncorrelated thus permitting any interaction 
type (predator-prey, mutualistic etc.) to exist in principle. 

We initiate the system with a set of individuals which are assigned random phenotype 
vectors. At each subsequent timestep a randomly selected individual is annihilated with 
probability Pkui = 0.2 whereby it is removed from the system and the single resource 
unit associated with it is returned to the resource bath, R{t). During the same timestep, 
another individual is randomly selected to reproduce with probability Prepro- This value 
is determined through the use of a weight function. 



H{T<^, t) = a, gT.rJ(T",T)n(T,0 _ y _ iV^_ 

'ET6rC(T",T)n(T,t) ^ ^ ) K ^) 3 ^^^^ K) 

This is monotonically mapped to the interval [0,1], appropriate for a probability mea- 
sure, by using the following function, 

_ exp[i/(T", t)] 

The sums of Eq.(2) are made over the points in phenotype space, T, and the occu- 
pancies (population associated with each phenotype), n(T,t) are used to account for the 
multiplicity of individuals with the same phenotype vector. We consider here a well mixed 
system of constant spatial size, although spatial extent is not explicitly considered. For 
clarity we reiterate at this point that the phenotype space is a pre-defined, complete set 
of all possible phenotypes and it is evolution and contingency that select the actualised 
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phenotypes in the evolved system. The J-matrix is similarly a pre-defined complete set 
of all possible interactions for all possible phenotypes that may exist in potentia. 

The correlation measure is used in the first term of Eq.(2) to restrict the impact of 
the interaction sum. It represents the fact that interactions are shared amongst members 
of the same species. For example, a tree may provide a volume of fruit to a solo member 
of a species but the provision must be shared with reduced efficacy if there are many 
members. So the overall effect of this denominator, on an individual, is to dampen its 
interaction sum as a whole thus representing the effect of distributing all interaction effects 
amongst the individual's own species members. Another example could be a wildebeest in 
proximity to a lion. The negative predatory effect (or predation probability) of the single 
lion on that specific wildebeest is decreased if there are many wildebeest about to select 
from. The interaction is damped by the presence of other members of the individuals own 
species. This aspect is not ubiquitous in species interactions but does feature in many 
cases. In recognition of the phenotypic variation inherent in a species, the sum over the 
correlation values, C(T", T^), accounts for the fact that species members have different 
but similar phenotypes. This can be seen as a species description in itself. 

When referring to species in this particular model we specifically mean the wild-type 
diversity in the sense that the mutant cloud is neglected in this categorisation. But, all 
interaction effects from all extant phenotypes, mutants included, are accounted for in 
every interaction sum. The mutant cloud itself is in fact very sparse as we have elected to 
use a low mutation rate, Pmut — 0.0002. As a consequence, the phenotype distribution is 
essentially a set of delta points of high occupancy with infrequent mutants existing with 
low occupancy. This makes the recognition of the wild-type diversity very simple as each 
species is massively dominated by the population of the wildtype. 

The second term of the weight function represents intra-specific competition and uses 
the same correlation measure as before. Similar but distinct phenotypes are likely to be 
in competition for resources, space etc. that are specific to their niche. The correlation 
measure accounts for this similarity. 

The third term represents competition for a conserved vital resource that all phe- 
notypes require for survival. Any successful reproduction event produces an offspring 
that assumes a unit of resource from the bath, R{t). The conservation requirement, 
R{t) + N(t) — constant, where N{t) is the global population of the system, means we 
have a carrying capacity for the system as a whole. It's functional form represents the 
number of system members competing per unit of available resource. 

The parameters, ai = 0.5, a2 = 0.01, 03 = 0.2 are the selection, conspecific competi- 
tion and resource competition parameters respectively. These are subjectively chosen to 
allow interaction controlled dynamics and a sufficient number of species to develop. A 
value of = 0.05 is used throughout the simulations. 

Upon successful reproduction, the individual produces a single offspring that assumes 
one unit of resource from the resource bath. As the bath depletes the probability of 
reproduction for any phenotype is reduced due to the third term in Eq.(2). The offspring 
phenotype is identical to that of the parent (clonal) unless, with probability Pmut — 
0.0002, a mutation occurs. This is effected by shifting a randomly selected trait value 
by a random amount V', that is gaussian distributed with fi = and a = ^, where 
^ = 250 is the correlation length of the phenotype space. Any mutated offspring will 
have an interaction set that is similar but not identical to its parent in accordance with 
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Figure 1: Degree distributions for the Tangled Nature model simulations. Shown are ensemble 
averaged data taken from all networks with number of species, D = {19,26,29} over 50 simu- 
lation runs of 10^ generations each. The exponential forms are highlighted by comparison with 
a binomial distribution of Z) = 29 and equivalent connectance, C~0.145 to the simulation data 
of the same number of species. 



the correlation. 



3 Results: degree distribution 

The degree distribution, P{k) represents the probability of a randomly selected species 
having k interactions with other members of the system. Given a network where edges are 
randomly distributed with identical probability, the degree distribution takes binomial 
form. Many ecological networks present distributions that are far from binomial, assum- 
ing forms more appropriately described as exponential or power law (Dunne, Williams, 
& Martinez, 2002)(Jordano, Bascompte, & Olesen, 2003). The version of the Tangled 
Nature model that incorporates correlated inheritance produces exponential distributions 
in nearly all cases (Fig.l). 

There are theoretical reasons why these forms, particularly the power law, should 
offer stabihty to such systems (Albert & Barabasi, 2002), so selection can be proposed as 
a causal effect for their observed existence. In a nutshell, the network structure permits 
a greater degree of resilience to random species extinctions. The stability arguments 
are certainly valid but it may actually be the case that the distributional forms appear 
as a consequence of the internal dynamics. When we perform simulations with random 
inheritance the degree distributions revert back to binomial form. This implies that 
the correlations are an essential requirement for constructing our exponential networks. 
To support this theory we shall now present a network model that produces a range of 
non-binomial distributions through correlated dynamics, without any form of selection. 
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4 Network evolution model 



We consider a system with a fixed number of species, D, eacli defined by generalised 
interactions with subsets of the other system members. Self-connections are excluded here 
and as the interaction types are not explicitly considered, the connections are regarded 
as undirected. This represents a simplified species interaction network that, in principle, 
embodies interaction types such as mutualism in addition to the usual food-web based 
predator-prey relationships. As a result, the networks we consider here are not expected to 
assume non-cyclic tree structures nor the stratified trophic levels associated with resource 
flow. 

We now initiate dynamics representing extinction and correlated speciation with the 
constraint of an invariant species number. Newly speciated members are seen to super- 
sede extinct members without any implication of cause nor effect. This invariance of 
the species number can be treated as a consequence of a carrying capacity and whilst 
simplistic is a reasonable approximation to an ecosystem. At a timestep, we randomly 
select with uniform probability a species and all its associated edges to be deleted from 
the network. A second parent species is then randomly selected from the remaining D—1 
and is duplicated in the form of a daughter. All species connected to the parent are now 
given connections to the daughter with probability Pg- All species unconnected to the 
parent are given connections to the daughter with probability, P„. The determination of 
an edge between the daughter and parent is made with probability Pp. 

At mean-field level the system admits a single attracting fixed point in the con- 
nectance, C^Cq, which is determined by the duplication probabilities (Eq.6). This can 
be shown when considering the following evolution equation for the number of network 
edges, 

Et+2 = Et-kt + Pekt+i + Pn{D-2-kt+i)+Pp, 
- IE - 2 - 

h = h^,^^—^{Et-h). (4) 

The second term accounts for the loss of edges when a node with mean degree, kt is 
deleted in stage one. The third and fourth terms then represent the increase in edges 
incurred at stage two when one of the remaining D — 1 nodes is duplicated. Imperfect 
duplication is made possible at this point by the inclusion of the probabilities, and 
Pn. The fifth term represents the contribution from the daughter-parent connection. 

As we only observe the system after every two steps we may rescale the time variable 
increment, (t + t')^(t + |). Thus Eq.(4) may be expressed as a recursion relation that 
has a general solution of the form, 

Et = Ail + aY-^. (5) 

a 

This admits a single fixed point dependent upon the parameter values which can be 
expressed in terms of the network connectancc, 

_ KjP - 2) + P, 

- D-l-(Pe-P.)p-2)- 

This point is an attractor for all physically relevant values; D>2; < Pe,Pn,Pp < 1- 
By using this relation we set the probabilities Pg, Pn dependently such that Pp=Co, so 
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allowing P^efCo, 1] to be used as a control parameter for the fidelity. As a result the 
dynamics range from fully random, Pe=P„=Pp=Co, generating binomial networks, to 
fully correlated, Pe=l, Pd=0, Pp=Cq, where daughter and parent are indistinguishable. 

4.1 Numerical results 

The ensemble averaged degree distributions for various fidelities are shown in Fig. (2). For 
the uncorrelated, random duplication process, Pd=Pn=Co we see binomially distributed 
networks as expected. As we increase the fidelity, the distributions exhibit longer tails 
and conform more to exponential than they do binomial. This progression continues 
until we have a process of perfect duplication at which point the distributions become 
power-law-like. 

Shown in Fig. (3), arc the ensemble averaged degree distributions for perfect dupli- 
cation networks of Co=0.01 at different system sizes. For small connectance values the 
distributions conform loosely to a power law. But as we reduce Co towards zero, the 
functional form becomes a discernible power law with exponent, 7~1. The networks 
achieved here are extremely sparse and the majority of the weight is held in the P{k=0) 
degree. So the form we see here represents (near) scale-free fluctuations. In the limit 
that Cq— a finite sized system asymptotically achieves the absorbing state of zero con- 
nectance. As we take the thermodynamic limit in the network size, though, we have a 
system where P{k=0) tends to unity and the remaining distribution achieves power-law 
form but with P{k>0) tending to zero. The simulations suggest that in these limits we 
have a network of zero connectance with 1/k scale-free fiuctuations in the degree. It can 
be verified analytically that this is indeed the case (Laird & Jensen, 2006). 

The power law distribution appears under the conditions of perfect correlation and 
taking the limit of vanishing connectance (or equivalently, Pp— >0). Whilst this is not 
representative of a real system we may relax these conditions whilst still achieving near 
power law forms. Due to the irresolute nature of field data, observed distributions are 
naturally indefinite anyway. If we assume though that they do conform to power-laws 
then it is instructive to compare their exponents. The low value, 7~1 is of significance 
here as ecological systems with power law distributions generally display such low range 
exponents, (Jordano et al., 2003) (Montoya & Sole, 2003). Biological systems in general 
exhibit exponents in the range, 7G(1, 2), (Chung, Lu, Dewey, & Galas, 2003) which can be 
reproduced with duplication-divergence growth models. Ecological degree distributions 
find themselves at the lower end of this range, as do those from our model, and the lack 
of network growth might be an important factor in achieving this. 

Field data degree distributions vary in form, with no obvious reason why they should 
take any one in particular, let alone why multiple forms do actually occur. Our model 
suggests that these distributions may simply arise as a consequence of correlated dynamics 
with the nature of the distribution being determined by the level of correlation. The 
idealised system modelled here fails to account for other defining processes, such as species 
invasion, and these will have an impact on the correlations. But, ultimately any species 
derives as an offshoot of an ancestral species so correlations will be inherent to the system 
even if a level of decorrelation is occurring. 
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Figure 2: Degree distributions for D = 200, Co = 0.01 systems produced using the imperfect 
duplication process. Prom short to long tail we have P(i={0.01, 0.25, 0.75, 0.95, 0.99, 0.999}. 
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Figure 3: Perfect duplication degree distributions for networks with number of species, 
L>={25, 50, 100, 200, 400, 800, 1600}, Co=0.01. Each has been normahsed with the exclusion 
of the k=0 support and rescaled to overlap for visual purposes. 



5 Connectance-species relationship 

We now return to the correlated Tangled Nature model and the observed inverse rela- 
tionship between the connectance, C and number of species, D. Connectance is a simple 
and commonly used measure of a network and it appears explicitly in the May-Wigner 
criterion of ecosystem stability (May, 1974). It represents the relative proportion of 
edges existing in a network when compared to the maximal set of all possible edges. 
Studies of field data have indicated that the network connectance and the number of 
species are inversely dependent. The networks created in the correlated version of the 
Tangled Nature model also exhibited a similar form of dependence (Fig. 4). For the real 
ecosystem, the constant connectance hypothesis and link species scaling law (Martinez, 
1992) (Cohen, Briand, & Newman, 1990) have been proposed to describe this relationship 
but neither viewpoints fully encompass the observed form. The May-Wigner criterion 
for stability naturally incorporates a, 1/D dependence for the connectance which is qual- 
itatively similar to the real system. Recently, an analytical population dynamics model 
that incorporates species invasion has provided an improved parametrised functional fit 
(McKane et al., 2000). The connectance value used in their work was that of the bare 
network which is analogous to the 6 parameter used in the Tangled Natiuc model (see 
section 2). By considering aspects of the Tangled Nature model we propose an alternative 
mechanism that recognises the role of connectance as a variable in the selection process 
and so account for the fact that it may deviate from the bare network value. 

We consider a hypothetical set of all viable organisms each described by distinct phe- 
notypical traits that determine their interactions with other phenotypes. Each organism 
interacts with a fixed subset of all other possible organisms which as a whole defines a 
complete species interaction network. An ecosystem represents a sub-network of this com- 
plete network where evolution has moulded the phenotype distribution such that a small 
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proportion of the possible organisms are extant. Realistically the strategies employed by 
organisms can depend upon population frequencies and environmental cues, so defining 
a complete superset of pairwise interactions is recognised to be an approximation. 

We may now ask the question: given a random set of biologically viable organisms 
what properties would the realised interaction network exhibit? If the complete set J has 
a connectance Cj with a binomially distributed degree then a purely random subset will, 
on average, have the same connectance and degree distributional form. For the following 
argument we shall assume that our complete bare network has such properties. 

In an evolving system, the effects of selection are partly determined by phenotypic 
interactions. Positive interactions are beneficial to an organism so we can conceive that 
selective processes incur an increase in such interactions. Dynamical systems theory sug- 
gests that mutualism should not be overly prevalent in community dynamics for stability 
reasons (May, 1974) but the hypothetical average interaction strength may still tend 
towards positive values. Mutualism clearly does occur in real systems and interactions 
described as parasitic may in reality be more commensal than antisymmetric (May, 1974). 
The networks exhibited in the Tangled Nature model attained higher than average con- 
nectance values with significantly greater mean interaction strengths. This demonstrated 
that selection processes applied at the level of individuals can cause a global increase in 
connectance. 

We can perceive this increase in edges as a shift upwards from the mean of the binomial 
distribution. But this shift is counterbalanced by the diminishing binomial probability 

of achieving higher edge numbers. The phase space of network configurations is heavily 
dominated by realisations that lie around the binomial mean so a selection-driven search 
in this space is unlikely to encounter high connectances. We estimate the shift as being a 
certain proportion, s, of the fluctuations one encounters in an unbiased binomial network, 
a{D, Cj). So in terms of the number of edges, E in the network we estimate. 



E = <Ej{D,Cj)>+sa{D,Cj) 
= ErnCj + s[E^Cj{l-Cj)]-2, 



(7) 



where E^ is the total number of edges in a maximally connected network and < Ej{D, Cj) > 
is the mean number of edges for networks of size, D with edge probability Cj. The pa- 
rameter, s is a measure of selection strength and is assumed to be independent of the 
system size. The connectance of the evolved network is defined as, C = which gives. 



C = Cj + s 
^ Cj + s 



Cj{l - Cj) 
Em 

2Cj{l-Cj) 
D{D-l) 



(8) 



Here we have used our definition of the maximal edge set, Em = \D{D — 1). Other 
definitions that incorporate self-interaction and directional edges could also be applied. 

So we have demonstrated that, with probabilistic restraint, a selective force driving 
the system to achieve greater numbers of interactions will incur an inverse dependence 
between connectance and the number of species. The assumption that selection incurs 
a net increase in interactions is justifiable, but the assumption of a binomial network 
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Figure 4: Plot of ensemble-averaged mean connectances, < C > against the number of species. 
Error bars represent the standard error. The lower dotted line marks the null system con- 
nectance, Cj = 0.05, which the evolved systems clearly surpass. The overlaid functional form 
is that given by Eq.(8) using the correct background connectance, Cj = 0.05 and with a value 
of, s = 5.5 for the selection parameter. 



structure is a significant limitation. The first part of this paper demonstrated how corre- 
lations between species can give rise to exponential-like, and power law distributions. The 
connect ance-species relationship might depend on the choice of background distribution 
so a fuller argument would require consideration of a more appropriate distribution. Even 
so, the underlying point remains valid. Selection is driving the system towards greater 
numbers of positive interactions and so higher edge numbers and equivalently higher con- 
nectances. But as we increase the number of nodes in our subnetwork the probability of 
achieving an elevated connectance diminishes, resulting in a decreasing functional form. 
If we compare Eq.(8) to the ensemble data acquired from the Tangled Nature model we 
see that the form is qualitatively appropriate. Fig. (4). With a background connectance of 
Cj — 0.05, the value used in the simulations, the fit is good but not ideal. The simulation 
networks presented exponential-like degree distributions though so a formulation based 
upon binomial networks could be responsible for this deviation at higher diversities. 

6 Discussion 

We have shown here that the exponential degree distributions of the correlated Tangled 
Nature model may be attributed to dynamical rather than selective processes. Our 
network evolution model dynamics generate distributions ranging from binomial through 
exponential to power-law which encompasses the Tangled Nature model results and many 
of the forms observed in real ecological systems. In the case of the power law distribution 
our exponent 7~1 compares well with the low values associated with ecological networks 
that take power law form. The network dynamics are appropriate given the type of system 
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but the model is idealised and ignores other determining factors, such as migration. The 
random introduction of species acts to decorrelate the system so future work would need 
to take account of such wider considerations. 

Several theories have been proposed to explain the inverse relationship between con- 
nectance and the number of species observed in field data ensembles. We have provided 
reasoning for this phenomenon that has no reliance on stability arguments nor allometric 
scaling due to resource flow. We argue that evolution pushes an ecology to acquire posi- 
tive interactions but is restrained due to the improbability of achieving high connectances. 
A C~l/D relationship arises as a consequence of this which is in close agreement with 
the simulation data taken from the correlated Tangled Nature model. The fit is not ideal 
and field data appears to have a larger exponent than our relation (Havens, 1992) (Pimm 
et al., 1991), but our formulation does rely on certain approximations. In particular the 
assumption of a binomial network is unrealistic. We have demonstrated that dynamics 
involving correlated inheritance have a strong influence on the form of the degree distri- 
butions. An evolutionary search process naturally flows into regions of the network space 
where the degree distributions are non-binomial and this reflects the inherent correla- 
tions. The phase space has non-trivial structure and this can influence the probability of 
locating higher connectance regions during the search process. Accordingly, the proba- 
bilistic reasoning has limitations, but the underlying point still has validity. Future work 
could attempt to account for the non-binomial networks that feature in ecology and so 
produce a more appropriate formulation. 

One may ask how our findings can be related to real biological systems. In principle 
this should be possible within chemostat-based microbial evolution experiments. Two 
microbial ecosystems with the same initial composition should, as a result of evolution, 
develop interaction networks with different connectances that are dependent on their 
species number. The correlations in the inheritance process should ensure they continue 
to exhibit similar long-tailed functional forms for their degree distributions. 
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